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LONG-TERM  GOALS 

The  Unites  States  Navy  is  evaluating  the  feasibility  of  constructing  and  deploying  very  large  mobile 
offshore  bases  (MOBs).  Such  platforms  may  be  1.5  km  long  and  300  m  wide.  These  structures  fall  far 
outside  the  scope  of  typical  maritime  design  experience.  Proper  design  depends  upon  the  ability  to 
characterize  ocean  waves  fields,  particular  nonlinear  features  like  wave  “long  crestedness”  and  wave 
“groupiness.”  A  long-term  goal  is  to  provide  such  characterizations. 

OBJECTIVES 

The  specific  objective  of  this  work  is  to  develop  and  apply  techniques  for  the  objective  estimation  of 
wave  long  crestedness  and  wave  groupiness  from  spaceborne  synthetic  aperture  radar  (SAR)  data  of 
ocean  surface  waves.  Spaceborne  SAR  imagery  offers  the  unique  combination  of  high  resolution 
(25  m)  and  large  swath  (100  km)  necessary  to  measure  large  scale  wave  coherence. 

APPROACH 

The  authors  listed  above  are  the  primary  contributors  to  this  work.  The  basic  approach  here  is  two 
fold.  First,  develop  techniques  to  objectively  measure  two  properties  of  wave  coherence:  ocean  wave 
crest  length  and  long  crestedness.  Second,  apply  these  techniques  to  spaceborne  SAR  imagery.  The 
difference  between  crest  lengths  and  wave  groupiness  observed  in  SAR  imagery  and  those  predicted  by 
linear  theory  is  a  measure  of  wave  coherence  caused  by  nonlinear  wave-wave  interaction. 

WORK  COMPLETED 

SAR  Wave  Measurement 

Ocean  waves  are  typically  measured  at  a  single  location  as  a  function  of  time.  However,  assessment  of 
the  spatial  coherence  of  waves  requires  measurement  of  the  two-dimensional  surface  wave  field. 
Spaceborne  SARs  offer  important  advantages  in  measuring  spatial  coherence:  (1)  They  have  sufficient 
resolution  (typically  25  m)  to  capture  ocean  wave  patterns.  (2)  Spaceborne  SAR  image  swaths  are 
typically  100  km  wide  or  larger,  providing  a  synoptic  view  of  the  wave  field.  (3)  The  global  coverage 
of  spaceborne  SARs  allows  the  acquisition  of  ocean  imagery  in  the  midst  of  storms. 
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The  study  of  mechanisms  for  SAR  wave  begins  in  1977  with  the  description  by  Elachi  and  Brown  [1], 
Their  approach,  with  a  considerable  refinement,  remains  valid  today.  Additional  work  by  many  others 
has  filled  in  the  gaps  [2-5].  According  to  Bragg  scattering  theory,  the  radar  cross  section  in  a  SAR  im¬ 
age  is  proportional  to  the  roughness  on  the  ocean  surface  at  the  radar  wavelength  projected  on  the  sur¬ 
face.  For  near-nadir  incidence  (20°  to  25°),  the  wave  imaging  is  dominated  by  “tilt  modulation”  for 
range  (cross-track)  traveling  waves.  Specifically,  the  wave-induced  modulation  of  radar  cross  section 
(image  intensity)  is  proportional  to  the  tilt  or  slope  changes  of  the  long  ocean  surface  waves.  There  is 
high  confidence  in  the  ability  of  a  SAR  to  faithfully  image  range-traveling  waves. 

Azimuth  wave  imaging  is  intrinsically  tied  to  how  a  SAR  image  is  formed  [6,  7].  A  SAR  relies  on 
Doppler  information  to  achieve  fine-scale  azimuth  resolution.  As  an  azimuth-traveling  wave  propa¬ 
gates  it  advects  the  surface  so  that  a  component  of  the  periodic  wave  orbital  velocity  is  moving  along 
the  SAR  line-of-sight.  This  causes  the  apparent  azimuth  positions  of  scatterers  to  be  displaced  in  the 
SAR  image.  The  periodic  ocean  wave  motion  moves  scatterers  in  such  a  way  as  to  concentrate  and 
dilute  the  apparent  density  of  the  scatterers  in  the  SAR  image  at  the  spatial  frequency  of  the  long 
waves.  This  renders  azimuth-traveling  waves  visible  in  SAR  imagery.  This  imaging  mechanism  is 
known  as  “velocity-bunching.”  The  velocity-bunching  effect  is  proportional  to  both  the  orbital  veloci¬ 
ties  of  the  long  waves  and  to  R/V  ,  where  R  is  range  from  the  SAR  platform  to  the  surface  and  V  is 
the  radar  ground  velocity.  When  waves  and  R/V  grow  large,  the  azimuth  shifts  can  become  large 
enough  to  smear  the  image  and  significantly  degrade  the  azimuth  resolution. 

A  useful  guide  to  the  magnitude  of  this  degradation  is  the  equation  Amin  =  ( R/V)H\ 12  where  /Lniin  is  the 

minimum  detectable  azimuth  wavelength  and  H  s  is  the  significant  wave  height  (SWH).  Clearly,  the 
higher  the  sea  state  and  the  larger  the  R/V  ratio  are,  the  larger  the  minimum  detectable  azimuth 
wavelength  is.  For  a  high  altitude  (800  km)  SAR  platform  like  ERS- 1/2  or  Radarsat,  R/V  equals 
120  s.  The  minimum  detectable  azimuth  wavelength  varies  from  170  m  to  380  m  for  SWHs  from  2  m 
to  an  extreme  of  10  m.  Thus,  even  under  extreme  conditions,  high-altitude  SARs  should  be  able  to  im¬ 
age  azimuth-traveling  waves  at  wavelengths  that  will  resonate  with  MOBs.  SAR  wave  imaging  at 
lower  altitudes,  such  as  the  200-km  shuttle  altitude,  offers  greater  fidelity  in  wave  measurement. 

Crest  Length  Measurement 

One  aspect  of  wave  coherence  is  wave  “long  crestedness.”  To  economically  assess  wave  long  crested¬ 
ness  requires  the  development  of  automated  procedures  that  produce  long-crestedness  parameters  from 
image  archives  with  minimal  human  intervention. 

The  challenge  of  arriving  at  a  statistical  description  of  long  wave  crests  is  non-trivial.  We  present  here 
a  rigorous  and  repeatable  procedure  for  identification  and  measurement  wave  crests  in  imagery.  This 
procedure  is  dependent  on  the  discretionary  selection  of  only  a  few  key  parameters.  Although  a  unam¬ 
biguous  formal  definition  of  long  crestedness  may  exceed  our  grasp,  it  is  possible  to  compare  observed 
long-crestedness  with  the  long-crestedness  expected  if  ocean  surface  waves  were  purely  linear.  Here, 
we  describe  this  procedure  and  apply  it  to  a  sample  SAR  image  from  Hurricane  Josephine  acquired 
during  the  Shuttle  Imaging  Radar-B  (SIR-B)  mission. 


The  extraction  of  crest  length  statistics  from  SAR  imagery  involves  two  basic  steps:  (1)  the  filtering  of 
a  SAR  image  so  that  crests  are  conspicuously  identified  as  separate  “blobs,”  and  (2)  the  measurement 
of  crest  length  from  these  blobs. 

Wavenumber  Domain  Filtering.  We  first  Fourier  transform  SAR  image  intensity  so  that  filters  are 
applied  in  the  wavenumber  domain.  Let  I(x,  y )  represent  SAR  image  intensity  as  a  function  of  range, 
x ,  and  azimuth,  y  ,  position.  We  represent  the  Fourier  transform  of  fix,  y)  by  F(kx,ky)  where 
kx  and  k  are  range  and  azimuth  wavenumber,  respectively.  If  we  let  km[n  represent  the  minimum 
wavenumber  at  which  we  want  information,  we  can  construct  a  high-pass  filter  //,  (. kx ,  ky )  such  that 
H i {k x , k y  ) =1  for  k  >  kmm  and  0  elsewhere,  where  k  is  wavenumber  magnitude.  Typically,  we  set 
kmm  =  2k  / 1600  r/m.  To  remove  the  remaining  noise  we  threshold  based  on  the  values  in  the  image 
spectrum,  F2(kx,ky )  .  We  define  a  threshold  level  by  T  =  /rs  +  Aas  where  fls  is  the  mean  spectral 
value,  <7S  is  the  standard  deviation  of  the  spectrum,  and  A  is  an  adjustable  parameter,  typically  set  to 
3.  This  second  spectral  filter  is  defined  as  H2(kx,k  )  =1  for  F1  (kx,ky)  >  Ts  and  0  elsewhere.  By  ap¬ 
plication  of  these  filters,  we  generate  a  new  Fourier  transform,  F\kx,k  )  . 

Image  Domain  Filtering.  We  produce  a  filtered  SAR  image  by  computing  the  inverse  Fourier  trans¬ 
form  of  F\kx,k  )  .  To  separate  crest  from  non-crest  areas,  we  applied  a  threshold  to  the  resulting  im¬ 
age.  If  jJLl  is  the  mean  of  the  filtered  image  and  <Jl  its  standard  deviation,  then  the  threshold  level  is 
7j  =  jltj  +  Bol  where  B  is  an  adjustable,  typical  set  to  1.  Portions  of  the  image  greater  than  threshold 
are  set  to  1,  otherwise  they  are  set  to  0.  Each  contiguous  area  of  l’s  is  a  “blob.” 

Blob  Identification  and  Measurement.  There  are  standard  image  processing  procedures  that  are  able 
to  find  contiguous  regions  or  blobs  in  an  image  and  label  them  sequentially.  We  then  fit  each  blob  to 
an  ellipse.  For  each  such  ellipse  we  determine  the  semi-major  and  semi-minor  axes,  the  eccentricity, 
and  the  orientation  of  the  ellipse.  Twice  the  semi-major  axis  is  the  crest  length.  Here  we  apply  the 
procedure  to  a  25.6 x  25.6  km  area  to  the  Hurricane  Josephine  SAR  imagery.  After  full  processing  of 
this  image,  we  obtain  a  list  of  all  the  wave  crests  in  the  image  and  their  dimensions. 

Wave-wave  interactions  conspire  to  bunch  waves  up  and  create  wave  crest  lengths,  which  differ  from 
what  might  be  expected  if  all  wave  frequency  components  acted  independently.  By  using  a  simulation 
technique,  we  are  able  to  compare  the  observed  crest  length  statistics,  with  similar  statistics  that  would 
be  generated  from  images  having  the  same  spectrum  but  whose  frequency  components  are  independent. 

We  begin  this  simulation  by  using  the  spectrum  computed  from  observed  image.  We  assume  that  this 
spectrum  is  but  one  realization  of  the  ensemble  mean  spectrum,  which  represents  the  surface.  In  this 
single  realization,  each  spectral  element  is  random  variable  having  a/2  -distribution  with  2  degrees  of 
freedom.  The  mean  of  this  distribution  is  the  ensemble  mean  spectral  value.  We  approximate  this  en¬ 
semble  mean  spectrum  by  smoothing  the  observed  spectrum.  The  aim  is  to  produce  new  simulated  im¬ 
ages  with  the  following  properties:  (1)  The  ensemble  of  mean  spectrum  of  similar  images  equals  the 
observed  image  spectrum.  (2)  Each  spectral  component  is  independent  of  the  others,  i.e.  the  phase  of 
each  spectral  element  is  independent  of  all  other  phases. 


Fig.  1:  The  number  of  wave  crests  encountered  longer  than  the  abscissa  per  square  kilometer. 

Results  indicate  that  very  long  wave  crests  are  more  common  on  the  ocean  than  would  be  expected  if 
wave  components  were  independent.  Fig.  1  shows  the  number  of  wave  crests  longer  than  the  abscissa 
that  will  be  found  on  average  in  one  square  kilometer.  The  higher  curve  was  computed  from  the  origi¬ 
nal  SAR  image,  while  the  lower  curves  were  the  result  of  simulations.  The  two  curves  represent  the 
observed  and  simulated  images.  From  Fig.  1,  we  may  infer  that  in  the  observed  image  one  is  likely  to 
encounter  two  wave  crests  longer  1000  m  in  a  10  km'  area.  In  the  simulated  images,  it  is  not  likely  that 
a  single  such  wave  crest  would  be  encountered. 

Wave  Groupiness  Measurement 

Wave  groupiness  is  the  likelihood  that  large  waves  will  group  together  in  a  series  [8].  Any  offshore 
platform  will  have  more  difficulty  if  especially  large  waves  are  grouped  together  possibly  setting  up  a 
resonance.  We  suspected  that  wave  groupiness  might  be  apparent  in  the  two-dimensional  autocorrela¬ 
tion  function.  Fig.  2A  represents  the  two-dimensional  autocorrelation  function  of  an  image  near  Hurri¬ 
cane  Josephine.  At  the  center,  or  zero  spatial  lag,  there  are  local  correlations  on  scales  of  a  few  hun¬ 
dred  meters.  This  corresponds  to  the  existence  of  ocean  waves  in  the  image.  Put  more  simply,  a  wave 
crest  at  any  particular  position  is  correlated  with  other  wave  crests  one,  two,  or  more  wavelengths 
away.  The  correlation  decreases  as  one  moves  farther  away. 

Upon  more  careful  examination  of  this  two-dimensional  autocorrelation,  groups  or  packets  of  waves 
are  visible.  This  is  a  consequence  of  wave  groupiness.  By  rotating  the  autocorrelation  function  shown 
in  Fig.  2A  so  that  the  crests  are  horizontal  and  collapsing  this  two-dimensional  representation  to  one 
dimension,  the  groups  become  more  apparent.  Fig.  2B  is  a  plot  of  the  one-dimensional  autocorrelation 
function.  Again  we  note  the  oscillations  near  zero  lag  associated  with  surface  waves.  The  envelope  of 
the  oscillations  is  the  signature  of  groupiness.  It  has  a  spatial  period  of  2.5  km.  Hence,  for  this  par¬ 
ticular  wave  field,  one  would  expect  especially  large  groups  of  waves  separated  by  2.5  km,  about  10 
wave  periods. 


Autocorrelations  computed  from  simulated  wave  images  based  on  linear  theory  did  not  exhibit  this 
groupiness  behavior.  This  is  additional  evidence  that  empirical  data,  rather  than  linear  wave  theory  is 
needed  to  estimate  wave  groupiness. 


Fig  2.  A  two-dimensional  autocorrelation  function  (A)  and  one-dimensional 

autocorrelation  function  (B). 


RESULTS 

Spacebome  SAR  images  offer  the  prospect  of  measuring  ocean  surface  waves  on  scales  and  at  resolu¬ 
tions  necessary  to  make  judgments  about  the  effect  of  waves  on  very  large  offshore  platforms.  This 
past  year  has  been  occupied  developing  the  repertoire  of  image  processing  tools  necessary  to  measure 
waves  coherence  properties  from  SAR  images.  Preliminary  evaluation  of  imagery  in  the  vicinity  of 
Hurricane  Josephine  indicates  that  waves  exhibit  both  longer  crest  lengths  and  more  groupiness  than 
would  be  expect  on  the  basis  of  linear  theory  alone.  It  is  clear  that  wave  imagery  will  be  necessary  to 
characterize  properties  of  wave  coherence.  At  JHU/APL,  we  have  acquired  SAR  ocean  wave  data  sets 
from  Hurricanes  Josephine,  Bonnie,  Danielle,  and  Edouard.  In  the  coming  year,  we  plan  to  compute 
the  crest  length  statistics  and  autocorrelation  lengths  for  these  data. 

IMPACT/APPLICATIONS 

The  procedures  presented  here  allow  SAR  wave  imagery  to  be  used  to  gather  statistics  on  wave  coher¬ 
ence.  These  in  turn  can  provide  realistic  wave  coherence  statistics  to  aid  in  the  design  of  large  off¬ 
shore  structures. 

RELATED  PROJECTS.  “Synthetic  Aperture  Radar  Imagery  of  the  Ocean  Surface  During  the 
Coastal  Mixing  and  Optics  Experiment,”  Donald  Thompson  and  David  Porter 
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